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Abstract 

We find all the perfect fiuid G2 diagonal cosmologies with the property that the 
quotient of the norms of the two orthogonal Killing vectors is constant along each 
fluid world-line. We find four different families depending each one on two or three 
arbitrary parameters which satisfy that the metric coefficients are not separable 
functions. Some physical properties of these solutions including energy conditions, 
kinematical quantities, Petrov type, the existence and nature of the singularities 
and whether they contain Friedman-Robertson- Walker cosmologies as particular 
cases are also included. 
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1 Introduction 



The study of exact solutions of Einstein's field equations for a perfect-fluid energy- 
momentum tensor is a difficult task even when we assume the existence of a two- 
dimensional abelian isometry group acting on the manifold. Indeed, for these spacetimes 
the Einstein equations constitute a system of highly non-linear coupled partial differen- 
tial equations with two independent variables which cannot be solved at present in the 
general case without making further hypotheses. 

As we are interested in studying models which could represent inhomogeneous cos- 
mologies we will assume that the isometry group two-dimensional orbits are spacelike 
everywhere (the so-called G2 on S2 metrics) and we will also make the usual hypothesis 
that the fluid velocity vector of the perfect fluid is irrotational. Under these hypotheses, 
the models can be classified into four classes depending on whether the orbits of the 
group are orthogonally transitive (which means that the set of two-planes orthogonal at 
each point to the group orbits are themselves surface forming) and on the existence of 
integrable Killing fields. The simplest (and in fact the most studied up to now) of these 
classes is the so-called diagonal cosmologies in which the group orbits are orthogonally 
transitive and both the Killing vectors are integrable (which in turn implies that they 
can be chosen mutually orthogonal everywhere) so that it can be seen that there exist 
coordinates {t, x, z} in which the metric take the Einstein-Rosen's form 

ds" = F{x, t) i-de + dx") + G(x, t)P(x, t)dy^ + ^^p^dz"" . (1) 
^ ^ P{x,t) 

The first of the known perfect-fluid solutions under these hypotheses not satisfying an 
equation of state for a stiff fluid was found by Wainwright & Goode some time 
ago (1980). This family depends on three arbitrary parameters and the perfect fluid 
satisfies a linear equation of state with the pressure proportional to the density. It 
has the particularity that the metric coefficient P depends only on t, or equivalently 
the three-spaces orthogonal to the fluid velocity are conformally flat. This solution has 
also the property that all the metric functions are separable in the coordinates t and 
X, that is to say, they are products of one function depending only on t and another 
function depending only on x. New metrics were later found in [Q], and 0, all of 
them being also separable metrics with explicit assumptions on the form of the functions 
involved. Finally, the case of separable diagonal cosmologies was extensively treated 
in a remarkable paper 0, where the general equations for separable cosmologies in 
co-moving coordinates (not including the family due to Wainwright & Goode nor the 
solutions satisfying the equation of state for a stiff fluid) were written and studied. 
Explicit families were also found showing that all kinds of singular behaviours were 
possible even when the energy conditions were satisfied everywhere (model with a big- 
bang-like singularity and/or big crunch singularity, metrics with timelike singularities, 
models with pressure and density regular everywhere but with singularities in the Weyl 
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tensor or complete models with no singularities at all). An extensive discussion of the 
relation between the existence of non-singular solutions satisfying energy and causality 
conditions and the powerful singularity theorems was published in |Q] . More recently, 
all the models with pressure equal to density and assuming separation of variables for 
diagonal cosmologies have been published in p. 

Thus, we see that the seek of exact solutions for diagonal cosmologies has been per- 
formed mainly assuming the Ansatz of separation of variables in co-moving coordinates. 
However, another approach has also been used in order to find new diagonal cosmologies. 
The use of the generalised Kerr-Schild transformation using Friedman- Robertson- Walker 



cosmologies as original base metric |]T0| has proven to be fruitful to find new diagonal so- 
lutions which are non-separable in the metric coefficients (some of these metrics have in 
fact only one Killing vector field). All of them, however, are of Petrov type D everywhere. 

Before stating clearly the assumptions we make in this paper, let us note that any 
Friedman-Robertson- Walker cosmology (which has a six dimensional group of isometrics 
acting on three-dimensional spacelike hypersurfaces) possesses also a two-dimensional 
isometry subgroup acting orthogonally-transitively on spacelike two-surfaces and with 
both Killing vectors integrable, so that the metric can be written in Einstein-Rosen's 
diagonal form. In fact, these metrics can be explicitly written as (see for instance []Tl| ) 



ds'^ = R\t) [-dr^ + dp^ + S2(p, k)d(t)'^ + ^'^(p, k)d. 



where R is an arbitrary function usually called scale factor, the function S(p, k) is given 
by 



sin p if k = 1 
S(p, k) = \ p if = 

sinh p if k = —1 



and the prime denotes derivative with respect to p. Thus, for the Friedman-Robertson- 
Walker cosmologies the function P in ([l|) is given by 

E'{p,ky 

which has the particularity that it does not depend on the coordinate r. Thus, the 
diagonal cosmologies such that the metric coefficient P in (|l]) does not depend on t 
obviously contain the Friedman-Robertson- Walker as a particular case and therefore it 
is worth considering such metrics without assuming other hypotheses on the form of 
the metric potentials. In fact, this is the aim of this paper where we find out all the 
diagonal G2 cosmologies with the property that the function P depends only on the 
spacelike coordinate. We find five families of exact solutions, each one depending on 
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arbitrary parameters. However, one of these families is in fact separable so that we do 
not study it here. All the other solutions are of Petrov type I and, consequently, they 



are not included in the non-separable Kerr-Schild families in |jlO|] so that they are, to 
our knowledge, new. For each family there exist ranges of the parameters such that the 
solutions satisfy energy conditions and also three of the four families contain Friedman- 
Robertson- Walker models for particular values of the arbitrary constants. One of the 
families represents a stiff fluid and the solution can be seen to be generated from a 



vacuum solution using the procedure due to Wainwright, Ince & Marshman 12 



The behaviour of the solutions is analysed and the existence of singularities and 
their nature is also studied. One of the families is complete and singularity-free, while 
other two have a big-bang-like singularity while they are regular in the future and in the 
spacelike directions. 

The plan of the paper is as follows. In section 2 we state the assumptions of this 
paper and give a coordinate independent characterisation, we also write the Einstein field 
equations and find the general solution of these equations. In section 3, we write down 
explicitly all the found solutions and analyse some of their physical properties, such as the 
energy density and pressure of the solution, whether they satisfy energy conditions, the 
Petrov type, the existence and nature of the singularities and the kinematical quantities 
of the fluid velocity vector. 



2 Field equations and their solutions. 

As stated in the introduction we will seek exact perfect-fiuid solutions of Einstein's 
field equations with an abelian two-dimensional group of isometrics acting orthogonally 
transitively on spacelike two-surfaces and such that both Killing vectors are integrable. 
Thus, there exist coordinates {t, x, y, z} in which the metric takes the Einstein-Rosen's 
form 

where the two Killing vector fields are given by 

dy' ^ dz 

Note that we could have set M = 1 and = 1 by means of a coordinate change of the 
type 

t' = t'{t), x' = x'{x), 

but in some cases is preferable to maintain these two functions in order to integrate some 
of the field equations and, in fact, this will be the case in our situation. 
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Our objective in this paper is to find out the co-moving perfect-fluid diagonal cos- 
mologies with the property that the function P of the metric does not depend on the 
variable t, that is to say, we will assume P{x). Rewriting this in a coordinate-independent 
way, we will flnd the perfect-fluid diagonal cosmologies which satisfy 

where u is the fluid velocity vector of the perfect fluid and and ^2 are two arbitrary 
mutually orthogonal Killing vector flelds. In order to see the equivalence of the two 
assumptions, we flrst note that in a perfect-fluid solution with a two-dimensional group 
of isometries acting orthogonally transitively on spacelike surfaces, the fluid velocity 
vector is necessarily hypersurface orthogonal and, in consequence, there exist coordinates 
{t, X, y, z} in which the metric takes the form with the fluid velocity vector given by 

u 

(co-moving coordinates). Except in the case that the function P in (0) is a constant 
(so that the metric possesses a three-dimensional group of isometries and is out of the 
interest of this paper), it is trivial to see that the only two mutually orthogonal Killing 
vector flelds are given by 

d , d 

oy oz 

where a and h are arbitrary non-vanishing constants. The condition (^ now reads 

^ d 





(p2) = ^ 9i (P) = 0, 



which is exactly the original assumption of P depending only on x. 

In order to write down the Einstein tensor (and all the other tensor components in 
this paper) we will use the orthonormal tetrad given by 



e'' = ^—dt, e^ = ^-dx, e^ = VGPdx, e^ = ^-dy. (4) 

so that the Einstein fleld equations for a co-moving perfect fluid energy-momentum 
tensor read (in units c = 8nG = 1) 

•S'oo = P, Su = S22 = 5*33 = p, Sal3 = {a ^ (3), 

where Sa/s (a = 0,1,2,3) stands fot the Einstein tensor and p and p are the energy 
density and the pressure of the fluid respectively. Obviously, there are combinations of 
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these equations involving the Einstein tensor alone and the calculation of the density 
and pressure can be performed once these equations are solved. In our case, the three 
equations which are non-identically satisfied are 

5*01 = 0, Sn — S22 = 0, 5*22 — 5*33 = 0, 

and they read explicitly 

1 1 / F P'^ P IF P P"^ 

5 _ L (^22 + 533) = — - -^M + tA'^M + -^M -tj^N + ^N+ 

^11 2 '^ 22 -r --^33; 2F V F F2 F G 2 F j? ^ j?2 ^ 

Li^n--^n'-^m--^m-^n + ^n--^n'-—n] =0, (7) 

F G 2F G 2G G G^ 2G J '^^ 

where the comma means partial derivative with respect to the variable indicated, the 
prime stands for ordinary derivative with respect to x and the dot indicates ordinary 
derivative with respect to t. 

In order to simplify this system of differential equations, it is convenient to define a 
new function K{x) through 

P' K 

— = -^. (8) 

We can assume K non-identically vanishing because we are not considering solutions 
with P constant. Thus, equation (H) can be integrated to give 

^^K + K' = O^Gix,t)- ^^^^ 



G ' ' ' K{xy 

where L{t) is an arbitrary no n- vanishing function depending only on the variable t. 
Substituting this expression for G into the equation we get 

FtK' F^L K'L 

In order to handle this equation some different cases have to be considered. First, let us 
assume L = so that we have 

Ft K' 
F K 
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As Ft = would imply that no function in the metric (Q) would depend on t and 
consequently dt would be a Killing vector, it follows necessarily that is a constant, 
but in this case the calculation of the density gives 

P = 7^ < 

and this solution is by no means physically reasonable. Thus, we must restrict our study 
to the case L ^ 0. If in addition we had K' = equation (P) would imply that the 
function F would not depend on x and then it can be easily seen that the metric 
possesses a third Killing vector field given by 



ox \ oz oy J 



so that we can discard this case. Thus we can assume K' ^ and L ^ and use the 
freedom that we still have in the coordinates t and x to fix these two functions locally to 

K = JsToe^ L = Koe' 

where we have introduced a constant Kq ^ because nothing prevents the function K 
to be negative. With this choice, equation (^) takes the simple form 

Ft - - F = 0, 

which can be trivially integrated to give 

F = e^H{t + x), (10) 

where H is an arbitrary non- vanishing function depending only on the variable u = t + x. 
Let us define a new function Z{u) by 

Z (11) 



H du 



so that (|^) takes the form 

(^ + ^)(M-iV) + |(M-iV')+^-ir2 = 0, (12) 

which is a linear first order ordinary differential equation for the function Z with the 
particularity (which in fact makes it difficult to study) that the coefficients do not depend 
on u but on other coordinates t and x. In order to handle this equation, we can consider 
Z as a function of the two variables t and x and impose the known condition that this 
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function depends only on the variable t + x. Thus, the equivalent system of partial 
differential equations for the function Z{t,x) is 



[Z^t + + z 



{M -N) 



Z,t - = 0. 



(13) 



In order to solve this system, let us first note that if M— vanished, we would necessarily 
have both M and constant and then the first equation in (|T^) would be clearly 
incompatible. Thus, we can define a new function ZQ{t,x) by the relation 



M-N 



(14) 



chosen in such a way that the system (|T3D takes the simpler form 

dtZo + d^Zo = e*+^- {2K^ - N') , 

M + N' 

r u X u ^ M — N 

We can evaluate the integrability condition dtd^Zo — d^dtZo = which gives 
{aK^ - N") {M-N) + (2K^ - N') (m + N') + 



(15) 



+Zoe-(*+^) In" + M + 



M-N 



0, (16) 



This equation suggests two possibilities, either with a vanishing or a nonvanishing coef- 
ficient of Zq. In the first case, the following two equations must be satisfied 

{aK^ - N") {M-N) + (2K^ - A^') (m + A^') = 0, 

(M - A^) (a^" + m) + N'"^ - = 0. (17) 

In order to find the general solution of this system of two ordinary differential equations 
that mixes functions of x and functions of t, we take the second derivative with respect 
to X and t of the second equation to get 

MN'" - MN' = 0, 

from what it follows that either M or N are constants or they satisfy 

N'" = bN', M = bM, 

where b is an arbitrary constant. Using this information, it is indeed very easy to find 
the general solution of (p!?!). There exist two different possibilities given by 
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Solution I 




No + Ky^ 



Solution II 

where A^^o ,-Ro and 5*0 are arbitrary constants. With these explicit forms for the functions 



M and A^, the partial differential equations ([TSD can be integrated to give Zq and therefore 



we can immediately calculate Z from its definition ( p!4D and H from (|TT]). The results 
are given by the following expressions (note that in the Solution I we have to distinguish 
between three subcases in order to perform explicitly the integral for H): 



Solution I with Sq = 0. 

^_ Rq-KI- ce-^(*+") 

Solution I with i?o = 0. 

Kge^(^+") + c 



H = Hq exp 



^ ^° {t + x) + ^e-2(*+^') 



2R 







H = Hq exp 



Solution I with i?o 7^ and Sq 0. 



2Sn - 2i?ne2(*+^) 



Solution II. 



-2ce-2(*+^) ^ H = Ho exp ice 



-2(t+x) 



where the arbitrary constants 6, c and Hq come from the integration of the differential 
equations (|T^ and (0). 



Thus, we have already found the general solution of the system (|1^) in the case that 
the integrability condition is identically satisfied. We can now move to the second case 
when we have 
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Now, the integrabihty condition (|T6|) gives exphcitly in terms of M and as 



iN" - 4^2) (M -N) + (N' - 2K^) (M + N') 

Zo = e*+^' (M - N) ^ , ' \ 1, (19) 

(M - A^) (A^" + Mj + N'^ - M2 

and substuting this expression in each of the equations in ([ISD we get two differential re- 
lations involving the functions M, A^ and their derivatives. The second of these equations 
reads explicitly 

(mm - M=^) A + (mm - mm) a + (mm - M^) A" + M (AA^' - A'A^) + 

M (A'A^' - A"A^) + M (A'A^" - AA^'") + AMK'^ {N'" - 2N") + 5 = 0, (20) 

where we have set 

A = N' -2K^, B = 2K^( N"'N' - 2NN'" - N"^ + 2N'N" + ANN" - AN'"^ 



The first equation in (|T5D gives an even more complicated relation between M and A^ 
but we will not need its explicit expression here. 

The study of the differential equation (|20|) is much more difficult that in the previous 
case and in fact some different subcases must be considered again. We will not make here 
the whole detailed study of these possibilities, and we will only indicate how to handle 
them in order to write down explicitly the solutions that arise and which also satisfy the 



first equation in ([T5|) . This is due to the fact that most of the calculations have been 
performed with the aid of computer algebra so that they are impossible to reproduce 
here. To begin with, we must note that if A vanishes identically, we have from ( plQ]) 
that Zq = and therefore we also have Z = and F = -FqC*, where Fq is an arbitrary 
positive constant. It is not difficult to see that in this case, the metric (0) possesses more 
than two Killing vector fields. Thus, we can assume that A is non-vanishing and divide 
equation (|20|) by A and take the derivative with respect to the variable x so that the 



term in (mM — M^) disappears. Consequently, the equation 

{MM - MM) (^)' + (MM - M^) )' + M J^^^l^J' + 



must be satisfied. From this equation we learn that either 

A' = bA, b = const 
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or we can divide the equation (pT]) by {A' /A)' and perform the derivative with respect 
to X so that the quadratic term (^MM — MM^ also disappears. It is now obvious that 
continuing with this procedure we finally get a linear third order differential equation 
for M 

aiM + + asM + a^M + as = 0, (22) 

where a^, (i=l,. . . ,5) are some complicated expressions depending only on x which con- 
tain some derivatives of quotients of the functions of x appearing in the original equation 
(pop. Except in the particular case that all the coefficients aj in (p2D are identically van- 
ishing it follows that M must satisfy a relation of the type 

M = aiM + azM + a^M + 04, (23) 

where 01,02,03 and 04 are now constants. This equation can be trivially solved for M 
although the explicit solution takes different forms depending on the multiplicity and 
the real or complex character of the roots of the polynomial 

X* - OiX^ - OaX^ - 03X - 04 = 0. 

Using these explicit forms for the function M some long but straightforward calculations 
allow us to restrict the values of 01,02, 03 and 04 which make the system (^) compatible 
and then find the solution for A^. Considering also the particular cases that we have 
been letting aside in the way to the equation (^) and which are not contained in ( p3D 
we finally get the following set of solutions which constitute the general solution of the 
system (|^) satisfying the condition (|ir 

Solution III 




-2t I c ^-6t 



5*16" 



Solution IV 



Solution V 




No + Sie 



where the constants Nq, Rq and 5*1 are arbitrary and the remaining constants 5*0 and c 
are restricted to be non-vanishing. 

Once the explicit forms the the metric functions M and are known, we can sub- 
stitute in ([ToD and back into ( |T4|) to find the expressions for the functions Z(t + x). 
Performing the integration of ([I^) we finally find the function H{t + x). The results are 
trivially found to be 
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Solution III Z = ^e2(*+-) ^ H = Ho exp ( ^e"^'+A . 

Solution IV Z= + ^h = Hod^'^^^ exp ( ^e'^'^A . 

Solution V Z = ^ ^H = i7oe^^*+"^ 

c + 4 " 

From the expression for H in the solution V we note that all the metric coefficients F, 
G and P are separable functions of the coordinates t and x and therefore this solution is 
contained in the general case of diagonal cosmologies separable in co-moving coordinates 
which was exhaustively studied by Ruiz&Senovilla @. Thus we will not study this 
solution here any further and we concentrate in the four families I, II, III and IV. 



3 Properties of the Solutions. 

In this section we will list all the solutions that we have found and give some of their 
physical properties such as the energy density and pressure of the perfect fluid, the 
kinematical quantities of the fluid velocity vector and the Petrov type. We will also 
discuss whether the solutions contain Friedman-Robertson- Walker models as particular 
subcases, if the energy conditions are satisfled for some ranges of the values of the 
parameters and the existence of singularities and their nature. 

For each of the four families of the solutions found in the previous section (Solutions 
I, II, III and IV) we have already written the functions M(t), N{x) anf H{t + x), while 
the metric coefficient G{x, t) is common for all the cases and is given by 

a = e'-\ 

In order to get all the metric components we still have to determine F{x,t) and P{x). 
The first of these functions is trivially obtained from H{x + t) using ([l0| ) and P{x) is 
found by solving the equation (|^) which reduces to a quadrature. 

For the general metric (0) with P depending only on x the non-vanishing components 
(in the tetrad (^)) of the kinematical quantities for the fluid velocity are given by the 
following expressions 



Expansion U 



Vf \F G 

Shear tensor au = 3-/^ ~ ~q ) ^22 = f^33 = "2^^^ 



1 Fr, 

Acceleration a = axdx = —dx 

2 F 
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while the vorticity obviously vanishes. The explicit form for the function G and the fact 
that F depends on the variable t + x except for a global factor e* imply the relation 



F ~ F G' 

In consequence, it is enough to list the acceleration and the expansion of the fluid for 
each solution in order to know all the kinematical quantities because the shear tensor 
components can be evaluated from 

1 iVm , , 

= cTsa = -2^11 = -3— (24) 



Before listing the possible non-separable families of solutions under the assumptions 
of this paper, let us note that some of the constants in M , and H above can be 
reabsorved by means of a linear coordinate change of the type 

X = aX + xq, 
t = aT + to- 

In the list that follows, this freedom has always been used to absorve all the spurious 
constants so that the parameters appearing in the solution are relevant (if ocasionally 
some irrelevant constants still remain, we explicitly state this fact as well as the reasons 
why this is so). 

Solution I 

As indicated in the previous section we have to distinguish between three subcases 
depending on whether some of the constants Rq or 5*0 vanish or not. 

Ro = 

The complete set of metric potentials for this solution is given by 

N = 1, M = 1 + ee-^^\ F = exp (at + (3a{t + x) + ec^ e'^''^'+''^) , 
G = e-^^*-^'), P = exp(2ce°^). 



where e is a sign. In this family, however, not all the constants a, (3 and c are true 
parameters identifying a particular solution. In fact, a change of the type (|25| ) allows us 
to fix a to any desired value whenever j3 is non-vanishing, while a is a true parameter 
when (3 = 0. However, we prefer to maintain the constants as they stand in order to 
study the cases /? = and P ^ together. 
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The energy density and pressure are given by 

ea^ (2/5 + 3) 



P = P 



4 exp (3at + (3a{t + x) + ec2e2«(*+^)) 



so that the perfect fluid satisfies a stiff equation of state. The energy (and pressure) is 
non-negative everywhere if and only if (3 and e satisfy 

e(2/3 + 3) > 0. 

The kinematical quantitites of the fluid velocity vector are 



e = ^- L a = - (2e^e^<'+^^ + (3) dx 

2 exp (ft + f a (t + x) + ef e2"(*+^)) 2 V 7 

The Petrov type is I except for the following particular cases 

• when c = and /? 7^ the solution is type D everywhere. 

• when c = 0, /? = the solution degenerates to conformally flat. Thus, the metric 
represents a Friedman- Robert son- Walker cosmology with density and pressure (we 
choose e = +1 so that the energy is positive) given by 

p = p = -a^e ^ . 
4 



50 = 

The list of metric coefficients is 

N=l + ee^"'-, M = 1, F = e"' exp (^^ (l - e(3^) {t + x) + ce-2'^(*+")) , 

Q ^ ^a(t-.)^ p _ I (e"^ + VT+^Y if e = +1 

I exp {(3 arcsin (e""^)) if e = — 1 

where e is a sign and the same considerations made in the previous case regarding the 
parameters a, [3 and c hold; when 3 — = the constants a and c identify uniquely a 
particular solution in the family while when 3 — e/?^ 7^ 0, the constant a can be fixed to 
any desired value by means of the coordinate change (^). 
The energy density and pressure read 
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V = P 



eo?c 



exp {^at +1(1- e/32) {t + x) + ce-2«(t+^)) 



so that again the perfect fluid obeys a stiff equation of state. The density (and conse- 
quently the pressure) is non-negative in the whole spacetime whenever we have 

ec > 0. 

The kinematical quantities of the fluid velocity vector take the form 

e = . ^ ^ ^, a = ^ (1 - e/?2 - 4ce-2»(*+-)) dx, 

4 exp [\at + f (1 - e/?2) {t + x) + |e-2«(<+^)) 4 V J 

and the Petrov type is I everywhere except for the following particular cases: 

• either when /9 = or when 1 — e/?^ = and c 7^ or when 9 — 6/3^ — and c = 
the metric is Petrov type D everywhere. 

• when 1 — el3'^ — and c = the solution is conformally flat. In this case the density 
vanishes and the solution is in fact the Minkowski spacetime. 

So # and ilo # 

The metric coeficients for this solution can be written in the following form 



Q^^a{t-x) p 



exp (c arcsin (e"^)) €1 — —1 



where both ei and 62 are signs and a can take the values +1, —1 or (obviously when 
any of the two signs ei or 62 are —1 then we necessarily have a — +1). The density and 
pressure are 

(rc2 (1 - U> - c,c') 

p = p- 



4 \eMt+x) _ 6162!'' exp (lat + ^ax - a (2b + eif ) {t + x)) ' 

which satisfy the energy conditions whenever the constants are restricted to 

€2 (4 - 46 - eic^) > 0. 
The kinematical quantities read 
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4 exp (fat + iax - a (6 + ei ^) (t + x)) |e2»(*+^) - eies ^ (e2«(*+^0 - ciCa) 



a e 



a 



,2a{t+x) 



[1 - eic^] + eac^ - 6162 + 46eie2 
4 (e2a{t+x) _ e^e^) 



and the Petrov type is again I except for the following particular values of the constants 

• when c = the solution degenerates to type D everywhere. 

• when eic^ = 1 and 6 = the solution is conformally flat and, therefore, it repre- 
sents a Friedman-Robertson- Walker cosmology with density (and pressure) given 
by (choosing 62 = +1 so that the energy conditions are fulfilled) 

p = p = a^e"^"*. 



Before finishing with the study of this Solution I let us indicate that all these p = p 
solutions can be seen to be generated from vacuum solutions using the procedure due to 



Wainwright, Ince & Marshman [T^. However, as far as we know, they were previously 
unknown. 

Solution II 

The metric coefficients are 

M = a + ee^"', N = a + e^"^", F = exp (at + ce-='"(*+")) , 



a. 



where e is a sign and the coordinate change (p5[) still allows us to fix the constant a 
equal to ±1 (depending on its sign) when a 7^ or set a equal to any desired value 
when a = 0. We allow this spurious constant in the family in order to consider the two 
subcases together. 

The expresions for the energy density and pressure read 

_ [^^^(Aat+2ax _ 4gce2'^* + 406^'^^) _ 2a2ee"* 

^ ~ 4 exp (3at + 2ax + ce-2«(*+^)) ' ^ ~ ^ ~ exp (ce-2»(*+^)) 
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so that the condition for positive energy density is, obviously 



3gg4at+2ax _ ^^^^2at ^ ^^^2ax y q_ ^25) 

In order to determine the ranges for the parameters a and c which make this condition 
possible everywhere we must distinguish between two cases depending on the sign of a. 
If a < it follows from the positivity of M that e = 1 and then the allowed ranges 
for the coordinates t and x are (we are choosing the constant a positive without loss of 
generality because a global change of sign in the coordinates t and x also changes the 
sign of a) 

i^>^log(-a), x>^log(-a). 
The condition ( P^D is fulfilled everywhere in this region whenever we have 

e = 1, c > , a < 0. 

4 

However, from the expression for the pressure it follows that there always exist regions 
of the spacetime which do not fulfil the condition 

p + p>0 

so that we will no longer consider this subcase as physically admissible. 

When a > 0, the coordinate x can take arbitrarily large negative values and the 
dominant term in the relation (^) when x —>■ — oo implies ec < 0. If e = +1 the 
coordinate t can also take arbitrarily large negative values and then ( pS]) in the limit 
t — > —oo clearly shows that the only possibility is c = 0. On the other hand, if e = — 1, 
the expression for M implies that a must be strictly positive so that we can use the 
freedom mentioned above to fix a = 1 and then the ranges for the coordinates t and x 
are given by 

t < 0, — oo < X < +00. 

It is not difficult to see that (|25|) is fulfilled in this region whenever c is restricted to 

3 

e = -l c>-. 

4 

The expression for the pressure shows that it is always bigger than the density and 
therefore positive everywhere. Thus, the only energy condition which is not fulfilled is 
p<p. 

It can be easily seen that when we approach t = along any timelike curve, the 
proper time remains finite. Consequently, either the spacetime is singular at t = or it 
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is extendible and our coordinate system only covers a portion of the whole inextendible 
spacetime. An analisys of the ciu^vaturc invariants near t = shows that all of them 
remain finite so that this solution is likely to be extendible. The same consideration holds 
when X ^ +00 where the proper length remains finite and all the curvature invariants 
have also a finite limit there. In order to show that this metric is in fact extendible, let 
us perform the change of variables 

t — — log [cosh (aT)] X — — log [sinh (ar)] 
Gj a 

where the range of variation of the new coordinates is given by 

T < 0, r > 0. 
The metric in this new coordinates take the following simple form 



cosh(2aT) 



,ccosh^(2aT)sinh^(2a.) ^ ^^2\ ^ ^^^^2 ^^2 ^ siuh (ar) 



where we have redefined a ^ 2a and we have renamed z ^ (}) and y ^ z because 
this metric is cylindrically simmetric with a regular axis of symmetry at r = where 
the so-called elementary flatness condition is satisfled. Thus, we need not extend the 
coordinate r as it represents a radial coordinate and the metric is already complete in 
the spacelike coordinates. Regarding the coordinate T, nothing prevents us to extend 
its variation range to the whole real line so that we have in fact extended the original 
spacetime which is included as the T < region. This solution is inextendible and the 
density and pressure in the new coordinates read 

a^ (4ccosh^(2aT) + 4ccosh2(2aT) sinh2(2ar) - s) 
P=- 



P = P + 



cosh(2ar) exp (ccosh^ (2ar) sinh^ (2ar)) 

8a2 

cosh(2ar) exp (ccosh^ (2ar) sinh^ (2ar)) ' 



which are now obviously non- negative (provided c > 3/4 as stated above). 
The kinematical quantities of this solution are 

asiuh(2rtr) (2r cosh'^(2(iT) f^mh^(2(ir) - s) 
^cosh(2aT) exp (f cosh^ (2ar) sinh^ (2ar)) ' 
a — accosh^(2aT) cosh(4ar)<ir. 

This solution is Petrov type 1 everywhere (provided c 7^ 0) except on the axis of symmetry 
r = where it degenerates to type D as requires the axial symmetry. 
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From the expressions above we see that both the density and pressure are regular 
everywhere and also the kinematical quantities are non-singular in the whole spacetime. 
In fact, it is easy to see that all the Riemann invariants of this solution are regular 
everywhere so that this metric contains no singularity at all. Thus we have found a cos- 
mological solution which satisfies the Strong Energy Condition everywhere and which 
contains no singularity at all. The first solution satisfying energy conditions and being 
complete and regular everywhere was found by Senovilla @], this solution was general- 
izad later to a rather large family in where diagonal separable metrics in co-moving 
coordinates were studied. Some recent work has tried to isolate that family as being in 



some sense unique |T^. However, very recently the author has found a solution which 
is also complete and non-singular (obviously satisfying energy conditions) and such that 
the line-element is non-diagonal and with separable metric functions |T^. The family 



that we present in this paper is the first known non-singular model with non-separable 
metric coefficients. 

Finally, let us note that in the particular case c = (which does not satisfy the energy 
conditions) all the kinematical quantitites of this solution except the expansion vanish 
identically and the perfect fluid satisfies a linear equation of state p = — 5/3p. Thus, 
this family contains a Friedman-Roberson- Walker model (with a physically unreasonably 
equation of state) when the constant c vanishes. 

Solution III 

The metric functions for this solution read 

M = 1 + ee-^"' + /5e-^"*, N = 1, F = exp (at + ec" e^^^'^'''^) , 
G' = e"(*-"), P = exp(2ce""), 



with energy density and pressure given by 

(See^'^* + 4ec2/3e2"(*+^) + 3/5) 2a^(3 
4exp(7at + ec2e2«(*+^)) ' ^ " ^ ^ exp (7at + ec2e2"(*+^)) ' 

From the form of the function M it follows that the coordinate t can take arbitrarily 
large values so that the expression for the density clearly imposes e = 1 in order to 
fulfil the energy conditions there. Similarly, the coordinate x can take arbitrarily large 
values and, given that the dominant term in the density when x — +oo is /^c^e^"'-*"'"^-', 
we have to impose /5 > in order to satisfy the non-negativity of the energy everywhere. 
Summing up, we restrict the constants of the solution to satisfy 

e = +l, /?>0 
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so that the density is positive everywhere. Unfortunately, the pressure is always bigger 
than the density except for /? = (this particular case is exactly the solution e = +1, 
/? = in the subcase i?o = of the Solution I, so that we will not consider it here any 
longer) . 

The ranges of variation for the coordinates t and x are 

— OO < t < +00, — OO < X < oo. 

The density and pressure diverge when t tends to — oo and it is very easy to show that 
this singularity is at finite proper time in the past. On the other hand, the density and 
pressure tend to zero when the time coordinate tends to +oo and it can be trivially seen 
that this happens at an infinite proper time in the future so that this solution presents 
a big-bang-like singularity in the past and it is complete in the future. It is also easy 
to prove that this solution contains no other singularities and that is is complete in the 
coordinate x. 

The kinematical quantities of the fluid velocity vector read 



a 



1 + e-2«* + /?e-6'^* J ^ ^, a = ac^e^'^^^+^^da;. 

4 ^ exp (|at + |e2'^(*+-)) 

It follows that the expansion is everywhere positive and it diverges when we approach 
the big-bang-like singularity while it tends to zero when t tends to +oo. It is interesting 
to note that the acceleration (and therefore also the sheartensor) and the Weyl tensor 
tend to zero near the big-bang and the density and pressure satisfy a linear equation of 
state in the limit near the big bang. In fact, the asymptotic behaviour of the density 
and pressure when t —>■ —oo is given by 

p ^a^/3e-^-, p ^«2^e-^- 

so that near the big bang they satisfy aproximately p = yp. Consequently, the metric 
is nearly homogeneous and isotropic (a nearly Friedman- Robert son- Walker geometry) in 
the vicinities of the big-bang and this geometry inhomogenizes and anisotropizes as the 
time goes on. In fact, the curvature singularity at t = — oo is Ricci dominated in the 
sense that 

hm ^-'^'^ , = 0, 

where obviously Cq,/3^5 stands for the Weyl tensor and Raji is the Ricci tensor. Fur- 
themore, in the particular case /? = (which implies the equation of state p = p) the 
singularity at t = — oo is an isotropic singularity in the sense that the metric is con- 



formally related to a metric which is regular at the t = —oo hypersurface (see [|T5| and 
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p6[] for a precise definition). Tlie Weyl tensor C^^^ in the natural coordinate cobasis 
{t,x,y,z} is regular but non-vanishing at t = — oo so that the so called Weyl tensor 
hypothesis is not satisfied. This hypothesis (see Penrose \T7\) assures that the appropri- 



ate thermodynamic boundary condition for the Universe is that the Weyl tensor should 
vanish at any initial singularity. It has been conjectured and partially proved by Tod 
(see |jl8|, |19[) that a perfect fluid model satisfying a barotropic equation of state with 
an isotropic singularity and satisfying the Weyl tensor hypothesis is necessarily an exact 
Friedman-Robertson- Walker cosmology everywhere. Thus, the particular solution (3 = 
here presented supports this conjecture due to Tod. 

Regarding the Petrov type of this solution, it can be easily seen that it is Petrov type 
I everywhere, except when c vanishes. In this particular case, the acceleration and the 
shear tensor also vanish and the density and pressure depend only on t so that they obey 
an equation of state. Thus, when c = 0, the solution is an exact Friedman- Rob er son- 
Walker geometry. 

Solution IV 

The metric coefficients of this solution are given by 

M=l + ee-^''\ N=l-e~^^^, F = exp (at + -a {t + x) + ec^e^''^'+'''> 



2 

2ace°'^dx 



e" 



-6ax 



where e is a sign and the expressions for the energy density and pressure are 
P=z 77^. — r. nrj—Tv^ P = P 



2 exp (fat + fax + ec2e2'^(*+^)) ' exp (^at + fax + ec2e2«(*+^)) ' 

The range of variation of the coordinate x is 

X > 

so that the positivity condition for the density when x +00 implies e = +1. With this 
value for e it follows that the p is positive everywhere and the pressure is also positive 
everywhere and satisfies 

p< p. 

Thus, this solution satisfies both the Strong and Dominant energy conditions everywhere. 
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Given that e — +1 the range of variation for the coordinate t is the whole real hne. 
When this coordinate tends to — oo, both the density and pressure diverge as 

This singularity at i — > — oo is at a finite proper time in the past and therefore this 
solution presents a big-bang- like singularity. On the other hand, when t aproaches -|-oo, 
the density and pressure tend to zero and the proper time also diverges so that this 
solution contains no singularity in the future. Regarding the coordinate x, the proper 
length diverges when this coordinate tends to +00 so that the metric is complete in that 
direction. However, the proper length of the apparent singularity at x = is finite so 
that either the solution is truly singular there or we must extend it beyond this value 
by means of an appropriate coordinate change. In order to perform this extension we 
define new coordinates X and T as 

i = 1 log [sinh(aT)] , 

X — — log [cosh(3aX)] , 
3a 

in which the metric takes the form 



ds'^ = sinh(ar) 



sinhi(ar) cosh5(3aX)e'='^^"^'("^)'^°^^^(^«^) [-dT^ + dX") 

^2ac f coshi(3aX)dX ^-Sac f cosh3(3aX)dX 

-dy H 1 az 



+ 



cosh3(3aX) cosh3(3aX) 

This metric is regular a,t X — (which corresponds tox — 0) and the range of variation of 
the new coordinate X can be extended to the whole real line. An analysis of this metric 
shows that the only singularity of this extended manifold is the big-bang singularity 
(which is now situated at T = 0) while the solution is complete in the future and in the 
spacelike coordinates. The density and pressure in the new coordinates read 

a^ (2c2 sinh^ (aT) cosht {3aX) + 9 sinh^ (aT) + 3 cosh^ (3aX)) 
^ 2 sinht (aT) cosht (3aX) e'=^ sinh2(aT) cosh§ (saX) 

4a2 

P = P- 



sinh2 (aT) cosh^ (3aX) e^'«i°h2(aT)cosh3(3aX) 



which are obviously positive everywhere and satisfying all the energy conditions at any 
point of the extended spacetime. Finally, the kinematical quantities of the fluid velocity 
vector are given by 
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^ a cosh(ar) (q + Ac'^ sinh^(ar) cosht (3aX)) 
4 sinh3 (aT) cosh^ (3aX)e^' sinh2(ar) coshi {3aX) ' 

^ a sinh(3aX) . sinh2(ar) coshi (3aX)) dX. 

4 cosh(3aX) V v y v 

This solution can be seen of Petrov type I independently of the values of the arbitrary 
parameters of the solution and, consequently, this family does not contain a Friedman- 
Robertson- Walker model as a particular case. 
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